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Abstract 

Optical conductivity of the weakly doped two-dimensional repulsive Hubbard model on the 
square lattice with nearest and next nearest hoppings is calculated within the generalized 
dynamical-mean field (DMFT+Sp) approach which includes correlation length scale ^ into the 
standard DMFT equations via the momentum dependent self-energy Sp, with full account of 
appropriate vertex corrections. This approach takes into consideration non-local dynamical cor- 
relations induced e.g. by short-ranged collective SDW-like antiferromagnetic spin fluctuations, 
which (at high enough temperatures) can be viewed as a quenched Gaussian random field with 
finite correlation length ^. The DMFT effective single impurity problem is solved by numerical 
renormalization group (NRG). We consider both the case of correlated metal with the bandwidth 
W and that of doped Mott insulator with [/ ^ ([/ — value of local Hubbard interaction). 
Optical conductivity calculated within DMFT-|-Sp demonstrates typical pseudogap behavior within 
the quasiparticle band in qualitative agreement with experiments in copper oxide superconductors. 
For large values of U pseudogap anomalies are effectively suppressed. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27+a, 71.30.+h, 74.72.-h 
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I. INTRODUCTION 



Pseudogap state is a major anomaly of electronic properties of of underdoped copper 
oxideaii^. We believe that the preferable "scenario" for its formation is most likely based on 
the model of strong scattering of electrons by short-ranged antiferromagnetic (AFM, SDW) 
spin fluctuations^. This scattering mainly transfers momenta of the order of Q = ^) (a 
— lattice constant of two dimensional lattice) leading to the formation of structures in the 
one-particle spectrum, which are precursors of the changes in the spectra due to long-range 
AFM order (period doubling) with non-Fermi liquid like behavior of spectral density in the 
vicinity of the so called "hot-spots" on the Fermi surface, appearing at intersections of the 
Fermi surface with antiferromagnetic Brillouin zone boundary (umklapp surface)^. 

In recent years a simplified model of the pseudogap state was studied^'^'^ under the as- 
sumption that the scattering by dynamic spin fluctuations can be reduced for high enough 
temperatures to a static Gaussian random field (quenched disorder) of pseudogap fluctua- 
tions. These fluctuations are defined by a characteristic scattering vectors of the order of Q, 
with distribution width determined by the inverse correlation length of short-range order 
K = and by appropriate energy scale A (typically of the order of crossover temperature 
T* to the pseudogap stated). 

It is also well known that undoped cuprates are antiferromagnetic Mott insulators with 
U ^ (f/ — value of local Hubbard interaction, W — bandwidth of non-interacting 
band), so that correlation effects are very important and underdoped (and probably also 
optimally doped) cuprates are actually typical strongly correlated metals. 

The cornerstone of the modern theory of strongly correlated systems is the dynamical 
mean- field theory (DMFT)^^^'^'^. At the same time, standard DMFT is not appropriate for 
the "antiferromagnetic" scenario of pseudogap formation in strongly correlated metals due 
to the basic approximation of the DMFT, which completely neglects non-local dynamical 
correlation effects. 

Different extensions of DMFT were proposed in recent years to cure this deficiency, such 
as extended DMFT (EDMFT)^°^^^, which locally includes coupling to non-local dynamical 
fluctuations, and, most importantly, different versions the so-called cluster mean-field theo- 
ries, such as the dynamical cluster approximation (DCA)^ and cellular DMFT (CDMFT)^. 
However, these approaches have certain drawbacks. First of all, the effective quantum single 
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impurity problem becomes rather complex. Thus, majority of computational tools available 
for the DMFT can be used only for small enough clusters^^, i.e. include mostly nearest- 
neighbor fluctuations. It is especially difficult to apply these methods to calculations of 
two-particle properties, e.g. optical conductivity. 

Recently we have proposed generalized DMFT+Sp approacbi^>i^ii^, which on the one 
hand retains the single-impurity description of the DMFT, with a proper account for local 
correlations and the possibility to use impurity solvers like NRGr^»2^, while on the other hand, 
includes non-local correlations on a non-perturbative model basis, which allows to control 
characteristic scales and also types of non-local fluctuations. This latter point allows for a 
systematical study of the influence of non-local fluctuations on the electronic properties and 
in particular provides valuable hints on physical origin and possible interpretation of results. 
Within this approach we have studied single - particle properties, such as pseudogap forma- 
tion in the density of states of the quasiparticle band both for correlated metal and doped 
Mott insulator, evolution of non-Fermi liquid like spectral density and ARPES spectrar^, 
"destruction" of Fermi surfaces and formation of Fermi " circs cis well as impurity scat- 
tering effects^. This formalism was also combined with modern LDA+DMFT calculations 
of electronic structure of "realistic" correlated systems to formulate LDA+DMFT+Ep ap- 
proach, which was applied for the description of pseudogap behavior in Bi2Ca2SrCu208^. 

In this paper we develop our DMFT+Sp approach for calculations of two-particle prop- 
erties, such as (dynamic) optical conductivity, which is conveniently calculated within the 
standard DMFT''--. We show that inclusion of non-local correlations (pseudogap fluctua- 
tions) with characteristic length scale ^ allows the description of pseudogap effects in longi- 
tudinal conductivity of the two-dimensional Hubbard plane. 

The paper is organized as follows: In section [TTl we present a short description of our 
DMFT+Ep approach. In section UTTl we derive basic DMFT+Ep expressions for dynamic (op- 
tical) conductivity, as well as formulate recurrence equations to calculate the p-dependent 
self-energy and appropriate vertex part, which take into account all the relevant Feynman 
diagrams of perturbation series over pseudogap fluctuations. Computational details and 
basic results for optical conductivity are given in section IIVI We also compare our results 
with that of the standard DMFT. The paper is ended with a summary section IVl including 
a short overview of related experimental results. 
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II. BASICS OF DMFT+Sp APPROACH 



As noted above the basic shortcoming of the traditional DMFT approacb^*^'^*^'^ is the 
neglect of momentum dependence of electron self-energy. To include non-local effects, while 
remaining within the usual "single impurity analogy" , we have proposecU^>i^ii^ the following 
(DMFT+Sp) approach. First of all, Matsubara "time" Fourier transformed single-particle 
Green function of the Hubbard model in obvious notations is written as: 

G{ie,p) = - , , N ^^T-T' £ = 7rT(2n + l), (1) 

where T,{ie) is the local contribution to self-energy, of DMFT type (surviving in the limit 
of spatial dimensionality d — > oo), while Sp(ie) is some momentum dependent part. This 
last contribution can be due either to electron interactions with some "additional" collective 
modes or order parameter fluctuations, or may be induced by similar non-local contributions 
within the Hubbard model itself. No doublecounting problem arises in this approach, as 
discussed in details in Ref.— . At the same time our procedure as stressed in Refs.— >^»^ 
does not represent any systematic 1/d expansion. Basic assumption here is the neglect of 
all interference processes of the local Hubbard interaction and non-local contributions owing 
to these additional scatterings (non-crossing approximation for appropriate diagrams)—, as 
illustrated by diagrams in Fig. [H 

The self-consistency equations of generalized DMFT+Sp approach are formulated as 
foUowaM^: 

1. Start with some initial guess of local self-energy T,{ie), e.g. Ti{ie) = 0. 

2. Construct Sp(ie) within some (approximate) scheme, taking into account interactions 
with collective modes or order parameter fluctuations which in general can depend on 
T,{iuj) and fi. 

3. Calculate the local Green function 

1 V 1 

Gu{^e) = ]^ E ,^ + ^_e(p)-S(2£)-Sp(z£)- 

4. Define the "Weiss field" 

g,\ze) = i:{zE) + GT,\te). (3) 
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5. Using some "impurity solver" calculate the single-particle Green function Gd{is) for the 
effective Anderson impurity problem, placed at lattice site and defined by effective 



action which is written, in obvious notations, as: 



Ses = - dri dT2Ci„{Ti)gQ\Ti - T2)cf^{T2) + drlJ n^{T)nii{T) . (4) 
Jo Jo Jo 

6. Define a new local self-energy 

sM = ^o~'M-Q'M- (5) 

7. Using this self-energy as "initial" one in step 1, continue the procedure until (and if) 
convergency is reached to obtain 

Gu{ze) = G4ie). (6) 

Eventually, we get the desired Green function in the form of ([1]), where T,{ie) and T,p{ie) 
are those appearing at the end of our iteration procedure. 



III. OPTICAL CONDUCTIVITY IN DMFT+Ep 
A. Basic expressions for optical conductivity 

To calculate dynamic conductivity we use the general expression relating it to retarded 
density-density correlation function x^i^y q}^'^: 

■ 2 

or(^) = -lim^x'^(^,q), (7) 

q^O 

where e is electronic charge. 

Consider full polarization loop graph in Matsubara representation shown in Fig. [2](a), 
which is conveniently (with explicit frequency summation) written as: 

<^{iu}, q) = ^ $i£i£'(i^, q) = XI ^) (8) 

ee' e 

and contains all possible interactions of our model, described by the full vertex part of 
Fig. El^b). Note that, we use slightly unusual definition of the vertex part to include the 
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loop contribution without vertex corrections, which shortens further diagrammatic expres- 
sions. Retarded density-density correlation function is determined by appropriate analytic 
continuation of this loop and can be written as: 

^ { [/(^+) - /(^-)] *f ^(q, ^) + ^(q, ^) - /(^+)*^^(q, ^) } , 

(9) 

where f{e) - Fermi distribution, e± = e±^, while two-particle loops ^f^{q,uj), $^^(q, cj), 
$^'^(q, uj) are determined by appropriate analytic continuations {ie + iu ^ e + u + i6, ie —>■ 
e ±i6,6 +0) in ([8]). Then we can conveniently write dynamic conductivity as: 

aH = lim (-^) J^Je{[f{e^) - /(._)] [«ff^(q,^) - $f^(0,^)] + 

[<l>f«(q, 0.) - <l>f «(0, - [$,^^(q, a;) - $^^(0, u)] } , (10) 

where the total contribution of additional terms with zero q can be shown (with the use of 
general Ward identities'^) to be zero. 

To calculate ^isi£'{iuj,q), entering the sum over Matsubara frequencies in (IHl), in 
DMFT+Sp approximation, which neglects interference between local Hubbard interac- 
tion and non-local contributions due to additional scatterings, e.g. by SDW pseudogap 
fluctuations^-, we can write down Bethe-Salpeter equation, shown diagrammatically in Fig. 
131 where we have introduced irreducible (local) vertex Uisie'{i(^) of DMFT and "rectangular" 
vertex, defined as in Fig. Et^b) and containing all interactions with fluctuations. Analytically 
this equation can be written as: 

e" 

where ^^^{iu, q) is the desired function calculated neglecting vertex corrections due to Hub- 
bard interaction (but taking into account all non-local interactions with fluctuations, con- 
sidered here to be static). Note that all g-dependence here is determined by ^^^{iLJ,q) as 
the vertex Uieie'iiuj) is local and g- independent. 

As is clear from (ITU]) , to calculate conductivity, we need only to find g^-contribution to 
^{iuj, q) defined in ([H]). This can be done in the following way. First of all, note that all the 
loops in fill I) contain g-dependence starting from terms of the order of g^. Then we can take 
an arbitrary loop (crossection) in the expansion of (fTT!) (see Fig. [3|), calculating it up to 
terms of the order of g^ , and make resummation of all contributions to the right and to the 
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left from this crossection (using the obvious left - right symmetry of diagram summation in 
Bethe - Salpeter equation), putting q = in all these graphs. This is equivalent to simple 
g^-differentiation of the expanded version of Eq. fllip . This procedure immediately leads to 
the following relation for g^-contribution to ([H]): 

^ hm ^ ^ ' ^ = 5^72^(z^,q = O)0°,(z^) 12 

e 

where 

__^tl(^^^q)_<(^ (13) 

with $°,(iu;, q) containing vertex corrections only due to non-local (pseudogap) fluctuations, 
while one-particle Green's functions in it are taken with self-energies due to both these 
fluctuations and local DMFT-like interaction, like in Eq. ([T]). The vertex 7j£(iu;,q = 0) is 
determined diagrammatically as shown in Fig. HI or analytically: 

7ie(zu;, q = 0) = 1 + ^ Uieie"{iuj)^ie"ie'{i^^, q = 0). (14) 



Now using Bethe-Salpeter equation (fTT!) we can write explicitly: 



-fie{lUJ, q = = 1 + > -g- = 

^ $^^(2^,q = 0) 

^ieie'{iuJ,q= 0) 



$o,(z^,q=0) • ^^^^ 
For q = we have the following Ward identity, which can be obtained by direct generalization 
of the proof given in Refs.— (see Appendix |A]) : 

{-iu)^is{iu, q = 0) = {-iuj) ^ ^ieie'iiuj, q = 0) = ^ G{ie + iu, p) - ^ G{ie, p). (16) 

e' p p 

Denominator of f[T^ contains vertex corrections only from non-local correlations (e.g. pseu- 
dogap fluctuations), while Green's functions here are "dressed" both by these correlations 
and local (DMFT) Hubbard interaction. Thus we may consider the loop entering the de- 
nominator as dressed by (pseudogap) fluctuations only, but with "bare" Green's functions: 

Go{ie,p) = — ^- — — — , (17) 

where T,{ie) is local contribution to self-energy from DMFT. For this problem we have the 
following Ward identity, similar to f|T6l) (see Appendix Rl): 

^ G{ie + iuj, p) G{ie, p) = $-e(icj, q = 0) [i:{ie + ioo) - S(ie) - iu] = 

p p 

= <^l{iu, q = 0) [AS(zcu) - H , (18) 



where we have introduced 

AE{iuj) = J:{ie + iuj) -Eiie). (19) 
Thus, using f|T6|) . ffTSl) in f|T5l) we get the final expression for •jieiiuj, q = 0): 



^i,{iuj,ci = 0) = 1 



Then ffT2l) reduces to: 



AS(icjl 



n 2 



(20) 



(21) 



Analytic continuation to real frequencies is obvious and using (]T2l) . ( 1211) in (iTOl) we can write 
the final expression for the real part of dynamic conductivity as: 
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Rea(ui) 



27r 



&[/(e_)-/(e+)]Re<'0f^(c. 



S^(e+) - S-^(5„) 
a; 

S^(5+) - S^(5_) 



(22) 



Thus we have achieved a great simplification of our problem. To calculate optical conductiv- 
ity in DMFT+Sp we only have to solve single-particle problem as described by DMFT+Sp 
procedure above to determine self-consistent values of local self-energies T,{e±), while non- 
trivial contribution of non-local correlations are to be included via (IT^ . which is to be 
calculated in some approximation, taking into account only interaction with non-local (e.g. 
pseudogap) fiuctuations, but using the "bare" Green's functions of the form (IT7I) . which 
include local self-energies already determined in the general DMFT+Sp procedure. Actu- 
ally fl22|) provides also an effective algorithm to calculate dynamic conductivity in standard 
DMFT (neglecting any non-local correlations), as (fT3!) is then easily calculated from a simple 
loop diagram, determined by two Green's functions and free scalar vertices. As usual, there 
is no need to calculate vertex corrections within DMFT itself, as was proven first considering 
the loop with vector vertices^"^. 



B. Recurrence relations for self energy and vertex parts. 

As we are mainly interested in the pseudogap state of copper oxides, we shall further 
concentrate on the effects of scattering of electrons from collective short-range SDW-like 
antiferromagnetic spin fluctuations. In a kind of a simplifled approach, valid only for high 
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enough temperatures^*^ we shall calculate T,p{iuj) for an electron moving in the quenched 
random field of (static) Gaussian spin fiuctuations with dominant scattering momentum 
transfers from the vicinity of some characteristic vector Q ("hot-spots" model^), using (as we 
have done in Refs.— i^^i^) slightly generalized version of the recurrence procedure proposed 
in^iii2i (see also Ref.=), which takes into account all Feynman diagrams describing the 
scattering of electrons by this random field. In general, neglect of fiuctuation dynamics, 
overestimates pseudogap effects. Referring the reader to earlier papers for details^'^'^'i^ii^, 
here we just start with the main recurrence relation, determining the self-energy: 

p) = A^-- . (23) 

Usually one takes the value of S^+i for large enough k equal to zero and doing recurrence 
backwards to k = 1 gets the desired physical self-energy S(ie, p) = T,i{ie, p)^'^>2i. 

In Eq. (1231) A characterizes the energy scale and k = is the inverse correlation length 
of short-range SDW fiuctuations, £fc(p) = ^(p + Q) and Vk = I'^p+ql + I'^p+ql for odd k 
while Skip) = £{p) and = |fp| + for even k. The velocity projections and 
are determined by usual momentum derivatives of the "bare" electronic energy dispersion 
e{p). Finally, s{k) represents a combinatorial factor, which here is always assumed to be 
that corresponding to the case of Heisenberg spin fiuctuations in "nearly antiferromagnetic 
Fermi-hquid" (spin-fermion (SF) model of Ref.-, SDW-type fiuctuations): 



As was stressed in Refs.— this procedure introduces an important length scale ^ not 
present in standard DMFT, which mimics the effect of short-range (SDW) correlations 
within fermionic "bath" surrounding the DMFT effective single Anderson impurity. 

An important aspect of the theory is that both parameters A and ^ can in principle be 
calculated from the microscopic model at hand^ -, but here we consider these as phenomeno- 
logical parameters of the theory (e.g. to be determined from experiments). 

Now to calculate optical conductivity we need the knowledge of the basic block $°g(ic<j, q), 
entering f|T3l) . or, more precisely, appropriate functions analytically continued to real frequen- 
cies: $°^"^(ct;,q) and ^^^^{u,q), which in turn define 0°^"^(u;) and (j)^^^{u) entering (l22l) . 




(24) 
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and defined by obvious relations similar to (fT3|) : 



= (26) 



By definition we have: 



$f^(.;,q) = 5^G^(^.+ ,p+)G-^(5_,p_)r^^(e„,p_;e+,p+) (27) 
p 

$f ^(c., q) = 5^ p+)G^(5-, P-)r^^(e-, P-; P+), 
p 

which are shown diagrammatically in Fig. [51 Here Green's functions G^{6^,p^) and 
G'^(£_,p_) are defined by analytic continuation (ze — > e ± i6) of Matsubara Green's func- 
tions (lH) determined by recurrence procedure (l23l) . while vertices r^^(£:_, p_; p_|_) and 
r^'^(e_, p_; £+, P4.), containing all vertex corrections due to pseudogap fluctuations are given 
by the recurrence procedure, derived first (for one-dimensional case) in Ref.~ (see also Ref.— ) 
and generalized for two-dimensional problem in Ref.— (see also Ref.-). The basic idea used 
here is that an arbitrary diagram for the vertex part can be obtained by an insertion of an 
"external field" line into the appropriate diagram for the self-energy2^>22i>2i. in our model we 
can limit ourselves only to diagrams with non-intersecting interaction lines with additional 
combinatorial factors s{k) in "initial" interaction vertices^'^'^. Thus, all diagrams for the 
vertex part are, in fact, generated by simple ladder diagrams with additional s(A;)-factors, 
associated with interaction lines^^i^ (see also^^). Then we obtain the system of recurrence 
relations for the vertex part T^^{e^, p_; £+, p+) shown by diagrams of Fig. [61 Analytically 
it has the following form^-, where we now also included contributions due to local (DMFT) 
self-energies, originating from DMFT+Sp loop: 

rffi(£_, p_; £+, p+) = 1 + Ah{k)Gt{e^, P-)Gf{e^, p+) x 

1 uj- £,(p+) + 5,(p_) - S«(5+) + S^(5_) - Sf:,,(e+, p+) + ^tii^-, P ) J 

xrf^(£_,p_;£+,p+), (28) 

and 

^(5±, p±) = _ , . , L.., ^ ^RA.- - ^ • (29) 



e± - £fc(p±) ± ikvkK - S^'^(£±) - S^;^(£±, p±) 
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The "physical" vertex T^^{e-, p_; e+, p+) is determined as T^^q^E-, p_; e+, p+). Recurrence 
procedure fl28|) takes into account all perturbation theory diagrams for the vertex part. For 
K ^ oo) (125]) reduces to the series studied in Ref.— (cf. also Ref.~), which can 
be summed exactly in analytic form. Standard "ladder" approximation corresponds in our 
scheme to the case of combinatorial factors s{k) in (125]) being equal to 1^. 

Recurrence procedure for r'^^(e_, p+; £+, p+) differs from (128]) only by obvious replace- 
ments A —>■ R and the whole expression in figure brackets in the r.h.s. of (128]) just replaced 
by 1: 

rf^i(£_, p_; £+, p+) = 1 + A''s{k)Gf{e^, P-)G^ie+, p+)r^^(£-, p; £+, p+). 

(30) 

Note that DMFT (Hubbard) interaction enters these equations only via local self-energies 
I]^'^(e±) calculated self-consistently according to our DMFT+Sp procedure. 

Equations (H]), ([23]), ([25]), ([30]) together with ([25]), ([26]) and ([22]) provide us with the com- 
plete self-consistent procedure to calculate optical conductivity of our model in DMFT+Sp 
approach. 

IV. RESULTS AND DISCUSSION 
A. Generalities 

In the following, we shall discuss our results for a standard one-band Hubbard model on 
a square lattice. The "bare" electronic dispersion in tight-binding approximation, with the 
account of nearest (t) and next nearest {t') neighbour hoppings, is given by: 

e{p) = —2t{cos p^a + cos Pyo) — it' cos p^a cos Pytt , (31) 

where a is the lattice constant. To be concrete, below we present results for t = 0.25eV 
(more or less typical for cuprates) and t'/t=-OA (which gives Fermi surface similar to those 
observed in many cuprates). 

For the square lattice the bare bandwidth is = 8t. To study strongly correlated 
metallic state obtained as doped Mott insulator we have used the value for the Hubbard 
interaction U = 40t and filling factors n = 1.0 (half-filling) and n = 0.8 (hole doping). For 
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correlated metal with W ^ U we have taken typical values like U = 4t, U = 6t and U = lOt 
for U > W. Calculations were performed for different fillings: half-filling {n = 1.0) and 
for n = 0.8,0.9 (hole doping). As typical values for A we have chosen A = t and A = 2t 
and for correlation length ^ = 2a and ^ = 10a (motivated mainly by experimental data for 
cuprates^i^) . 

To solve an effective Anderson impurity problem of DMFT we applied a reliable numeri- 
cally exact method of numerical renormalization group (NRG)^i2^, which, actually, allowed 
us to work with real frequencies from the very beginning, overcoming possible difficulties 
of performing analytical continuation numerically. Calculations were performed for two 
different temperatures: T = 0.088t and T = 0.356t. 

All necessary integrations were done directly, e.g. over the whole Brillouin zone (with the 
account of obvious symmetries), or wide enough frequency range. Integration momenta are 
made dimensionless in a natural way with the help of the lattice constant a. Conductivity 
is measured in units of the universal conductivity in two-dimensions: o'o = ^ = 2.5 10~^ 
Ohm~^. 

B. Optical conductivity in standard DMFT 

Optical conductivity was calculated for different combinations of parameters of the model. 
Below we present only a fraction of our results, which are, probably, most relevant for copper 
oxides. We shall start with presenting some typical results, obtained within our formalism 
in conventional DMFT approximation, neglecting pseudogap fluctuations, just to introduce 
the basic physical picture and demonstrate the effectiveness of our approach. 

Characteristic feature of the strongly correlated metallic state is the coexistence of lower 
and upper Hubbard bands splitted by the value of ~ t/ with a quasiparticle peak at the 
Fermi level^"^. For the case of strongly correlated metal with W > U we observe almost no 
contribution from excitations to upper Hubbard model in optical conductivity, as can be seen 
in Fig. [7] (where we show the real part of conductivity Iiea{uj)). This contribution is almost 
completely masked by a typical Drude-like frequency behavior, with only slightly non- 
monotonous behavior for u; ~ f/, which completely disappears as we rise the temperature. 

Situation is different in doped Mott insulator with U ^ W. In Fig. [H]we clearly observe 
an additional maximum of optical absorption for ~ f/, however, at smaller frequencies we 
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again observe typical Drude-like behavior, slightly non-monotonous for small frequencies 
due to quasiparticle band formation (see insert in Fig. [8]). 

These and similar results are more or less well known from the previous studies^*^, and 
are quoted here only to demonstrate the consistency of our formalism and to prepare the 
reader for new results, showing pseudogap behavior. 

C. Optical conductivity in DMFT+Sp 

1. Correlated metal 

Let us start the discussion of results obtained within our generalized DMFT+Sp approach 
for the case of > f/. 

In Fig. [H] we show our DMFT+ Sp results for the real part of optical conductivity for 
correlated metal {U = At) for two values of temperature, compared with similar data without 
pseudogap fluctuations (pure DMFT). We clearly observe formation of typical pseudogap 
(absorption) anomaly on the "shoulder" of Drude-like peak, which is partially "filled" with 
the growth of temperature. This behavior is quite similar to "mid-infrared feature" that 
is observed in optical conductivity of cuprate superconductors2Ii^. In Fig. [10] we show 
the behavior of Rea{uj) for different values of the pseudogap amplitude A. We see that 
pseudogap anomaly naturally grows with the growth of A. Fig. [TT] illustrates the dependence 
of Recr(u;) on correlations length of pseudogap (AFM, SDW) fluctuations. Again we observe 
the natural behavior — pseudogap anomaly is "filled" for shorter correlation lengths, i.e. 
as fluctuations become more short-ranged. At last, in Fig. [12] we demonstrate dependence 
of pseudogap anomaly in optical conductivity on correlation strength, i.e. on Hubbard 
interaction U. It is seen that the frequency range, where pseudogap anomaly is observed 
becomes narrower as correlation strength grows. This correlates with general narrowing of 
the pseudogap anomaly and spectral densities with the growth of correlations, observed in 
our previous work— li^. For large values of U pseudogap anomaly is practically suppressed. 
This is the main qualitative difference of the results of the present approach, compared to 
our earlier work— on optical conductivity in the pseudogap state. Comparing the data of 
present work for U = with similar data of Ref.— , it should be noted, that in this earlier 
work we have performed calculations of dynamic conductivity only for T = and used 
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simplified expressions, neglecting RR, AA-loops contributions to conductivity, as well as 
small frequency expansion^^, just to speed up calculations. These simplifications lead to 
some quantitative differences with the results of present work, where all calculations are 
done exactly using the general expression fl22]) . though qualitatively the frequency behavior 
of conductivity is the same. 

2. Doped Mott insulator 

Now we shall discuss our results for the case of doped Mott insulator with U ^ W . This 
case has no direct relevance to copper oxides, but is interesting from the general point of 
view and we present some of our results. 

The real part of optical conductivity for the case oiU = 40t is shown in Figs. [TSfT^ 

In Fig. [12] we show Re(j(c(j) for doped Mott insulator in DMFT+Sp approach several val- 
ues of pseudogap amplitude A. Obviously enough, pseudogap fluctuations lead to significant 
changes of optical conductivity only for relatively small frequencies of the order of A, while 
for high frequencies (e.g. of the order of [/, where the upper Hubbard band contributes) we 
do not observe pseudogap effects (see insert in Fig. [T3l) . For small frequencies we observe 
pseudogap suppression of Drude-like peak, with only a shallow anomaly for ~ A, which 
just disappears for smaller values of A or shorter correlation lengths. 

In Fig. [14] we show the similar data for the special case of t' = and n = 1, i.e. 
at half-filling (Mott insulator) for different values of inverse correlation length k = 
Conductivity at small frequencies is determined only by thermal excitations and pseudogap 
fluctuations suppress it significantly. Shorter correlation lengths obviously lead to larger 
values of conductivity at small frequencies. Transitions to upper Hubbard band are not 
affected by these fluctuations at all. 

V. CONCLUSION 

The present work is the direct continuation of our previous work— ^i^ii^, where we have 
proposed a generalized DMFT+Sp approach, which is meant to take into account the impor- 
tant effects of non-local correlations (in principle of any type) in addition to the (essentially 
exact) treatment of local dynamical correlations by DMFT. Here we used a generalized 
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DMFT+Sp approach to calculate dynamic (optical) conductivity of two-dimensional Hub- 
bard model with pseudogap fluctuations. Our results demonstrate, that pseudogap anoma- 
lies observed in optical conductivity of copper oxides can, in principle, be explained by this 
model. Main advantage in comparison to the previous work^^ is our ability now to study the 
role of strong electronic correlations, which are decisive in the formation of electronic struc- 
ture of systems like copper oxides. In fact, we have demonstrated an important suppression 
of pseudogap anomaly in optical conductivity with the growth of correlation strength. 

As we already noted in Ref.— qualitatively similar results on pseudogap formation in 
single-particle characteristics for the 2d Hubbard model were also obtained within cluster 
extensions of DMFT— li^. However, these methods have generic restrictions concerning the 
size of the cluster, and up to now were not widely applied to calculations of two-particle 
properties, such as general response functions, and, in particular, to calculations of dynamic 
(optical) conductivity. 

Our approach is free of these limitations, though for the price of introduction of additional 
(semi) phenomenological parameters (correlation length ^, and pseudogap amplitude A). It 
is mauch less time consuming, thus its advantage for calculations of two-particle response 
functions is obvious. It also opens the possibility of systematic comparison of different 
types of non-local fluctuations and their effects on electronic properties, providing more 
intuitive way to analyze experiments or theoretical data obtained within more advanced 
schemes. Note, again, that in principle both ^ and A can be calculated from the original 
model^. Our scheme works for any Coulomb interaction strength U, pseudogap strength 
A, correlation length ^, filling n and bare electron dispersion e(k). 

The present formalism can be easily generalized in the framework of our recently pro- 
posed LDA+DMFT+Sp approach, which will allow to perform calculations of pseudogap 
anomalies of optical conductivity for realistic models. It can also be easily generalized to 
orbital degrees of freedom, phonons, impurities, etc. 
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APPENDIX A: WARD IDENTITIES 

In this Appendix we present derivation of Ward identities used in the main text. Let us 
start with the general expression for variation of electron self-energy due to an arbitrary 
variation of the complete Green's function, which is valid for any interacting Fermi systen>2^: 

ASp = 5^f/pp,(g)AG'p., (Al) 
p' 

where Upp'{q) is an irreducible vertex in particle-hole channel, and we use 4-dimensional 
notations p = {ie, p), q = {i(jJ, q) etc. In the following we take: 

AEp = E+ - S_ = p+) - E{ie^, p_) (A2) 

and (in the same notations): 

AGp = G+-G^ = (G+G_)p(ASp - A{Go\), (A3) 

where A{Gq^)p = Gq^ — GqI, and the last expression was obtained using the standard Dyson 
equation. 

Note the similarity of Eq. (lAll) to the Ward identity for non-interacting electrons in the 
impure system, derived in Ref.— . 

Now substituting the last expression in flA3l) we get: 

ASp = J2Upp'{q){G+G^)piAEp, - A{Go%,). (A4) 
p' 

Iterating this equation we obtain: 

AEp = J2UpAG+G^)A-A{G,%>) + 
p' 

+ J2Upp"{G+GAp"Up"p'{G+GAp'{-HGo\') + ... . (A5) 
pi'pi 
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Multiplying both sides of (lASP by (G'+G'_)p and adding 

5^(G+G_),V(-A(Go^),0 = {G^G^U-A{G^%) 
p' 

we have: 

p' 

+{G+G-)p Upp//{G+G-)p'/Up'/pi{G+G-)pi + ...](— A(Gq = 
p" 

= Y.%p'iq){-A{G,%), (A6) 

p 

where is the complete two - particle Green's function determined by the following 

Bethe - Salpeter equatioii^^: 

%p'iq) = iG+G^)p6pp> + {G+G^)pJ2Upp'^pp'iq)- (A7) 

p' 

Finally we obtain: 

AGp = J2%P'm-MGo%), (A8) 
p' 

which the general form of our Ward identity. 

Summing both sides of ( lASj) over p and taking q = we obtain the identity (fT6ll used 
above. Similarly, taking the "bare" Green's function ( ITTI) we obtain (ITSll . 
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FIGURES 




FIG. 1: Typical "skeleton" diagrams for the self-energy in the DMFT+Ep approach. The first 
two terms are examples of DMFT self-energy diagrams; the middle two diagrams show some 
contributions to the non-local part of the self-energy (e.g. from spin fluctuations) represented as 
dashed lines; the last diagram (b) is an example of neglected diagram's leading to interference 
between the local and non-local parts. 




FIG. 2: Full polarization loop (a) with vertex part, which includes free-electron contribution in 
addition to the standard vertex, containing all interactions (b) . Here p± = pib3,£-|- = £zb^. 
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FIG. 3: Bethe-Salpeter equation for polarization loop in DMFT+Sp approach. Circle represents 
irreducible vertex part of DMFT, which contains only local interactions, surviving in the limit of 
d — > oo. Unshaded rectangular vertex represents non-local interactions, e.g. with SDW (pseudo- 
gap) fluctuations, which is defined similarly to Fig. EJb). 




FIG. 4: Effective vertex 7j£(ia;,q = 0) used in calculations of conductivity. 




FIG. 5: Diagrammatic representation of $°^"^(6<;,q). 
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FIG. 7: Real part of optical conductivity for correlated metal (U = 4t, t' = —OAt, t = 0.25 eV) in 
DMFT approximation for two values of filling factor: n = 1 and n = 0.8. Temperature T = 0.088t. 
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FIG. 8: Real part of optical conductivity for doped Mott insulator {U = AOt, t' = -OAt, t = 0.25 
eV) in DMFT approximation. Filling factors are: n = 0.8 and n = 0.9, temperature T = 0.088t. 
Small frequency behavior is shown in more details at the insert. 
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FIG. 9: (Color online) Real part of optical conductivity for correlated metal (U = At, t' = —QAt, 
t = 0.25 eV) in DMFT+ Sp approximation for two different temperatures: T = 0.088t and 
T = 0.356t. Pseudogap amplitude A = t, correlation length ^ = 10a, filling factor n = 0.8 
electrons per atom. 
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FIG. 10: (Color online) Real part of optical conductivity for correlated metal {U = 4t, t' = —OAt, 
t = 0.25 eV) in DMFT+Sp approximation — A dependence. Parameters are the same as in Fig. 
[21 but data are for different values of A = 0, A = t, A = 2t, and temperature T = 0.088t. 




FIG. 11: (Color online) Real part of optical conductivity for correlated metal {U = 4t, t' = —OAt, 
t = 0.25 eV) in DMFT+Sp approximation — dependence on correlation length. Parameters are 
the same as in Fig. El but data are for different values of inverse correlation length k = 
Ka = 0.1 and ku = 0.5, and temperature T = 0.088t. 




FIG. 12: (Color online) Real part of optical conductivity for correlated metal in DMFT+Sp 
approximation — U dependence. Parameters are the same as in Fig. [UJ but data are for different 
values U: U = 0, U = 4t, U = 6t, U = Wt and U = 40t. Temperature T = 0.088t. 




FIG. 13: (Color online) Real part of optical conductivity for doped Mott insulator (U = 40t, 
t' = -OAt, t = 0.25 eV) in DMFT+Sp approximation for different values of A = 0, A = 
A = 2t, and temperature T = 0.088t. Correlation length ^ = 10a, filling factor n = 0.8. Insert: 
conductivity in a wide frequency interval, including transitions to the upper Hubbard band. 
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FIG. 14: (Color online) Real part of optical conductivity for doped Mott insulator (U = 40t, 
t = 0.25 eV, t' = 0) in DMFT+Sp approximation for different values of inverse correlation length 
K = Ka = 0.1 and na = 0.5, temperature T = 0.356t and filling n = 1. 



